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brane Electromotive Forces and of Flows of Electrolytic 
Solutes 
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An analysis is given for the treatment of membrane transport phenomena in accord 
with the theory of steady state thermodynamics. A linear macroscopic theory for discon- 
tinuous systems is applied as a postulate. It is shown that as a consequence of the trans- 
formation properties of the Onsager reciprocal relations the definition of a membrane electro- 
motive force gives corollary definitions of the flows of electrolytic solutes as a whole in the 
form of hnear combinations of the flows of ionic constituents. It is shown that established 
conventions which set the activity coefficients of ionic constituents equal to unity at the 
reference state of infinite dilution lead to a particular definition of the membrane electro- 
motive force which may be apphed at any concentration. 



1. Introduction 

This paper reports the results of a part of a study 
directed toward the detailed apphcation of the 
thermodynamics of steady state processes to the 
investigation of transport phenomena at junctions 
between electrolytic solutions. The work was un- 
dertaken as a part of a research project under the 
sponsorship of the Office of Saline Water of the U.S. 
Department of the Interior. This project had the 
aim of improving the methods of measuring and 
reporting the electrochemical characteristics of 
membranes. 

In this paper we are concerned with the problem 
of treating permeability characteristics during steady 
states involving the transfer of electric charge without 
introducing the classical uncertainty involving the 
electrostatic potential difference between two chemi- 
cal phases of different composition [1, 2, 3, 4].^ The 
existence of such an electrostatic potential difference 
is postulated in discussions of the application of 
steady state thermodynamics to electrolytic trans- 
port problems. The usual practice [5] is to introduce 
the definition of the differential of the chemical 
potential of an ionic constituent in the form 



d^it=(dfit)a^=o+^{^d^ 



(1) 



where dfXi is the total differential of the chemical 
potential; z{ is the electrovalence ; J^is the faraday; 
and d^ is the differential of the electrostatic potential. 
The first term on the right, the ^'nonelectrical part/' 
may be put formally 

(dfjLf)ai^o=-Si'dT+VtdP+RTd[nmt+RTd[nYi 

(2) 

where >S1, Vi, m% and yt are, respectively, the partial 
molar entropy, the partial molar volume, the 



1 Figures in brackets indicate the literature references at the end of this paper. 



molality and the activity coefficient at constant 
temperature and pressure of the ionic constiuent. 
We note here that if ion constituents are defined as 
Maclnnis [6] defines them, mi and hence din m\ are 
operational quantities without extra-thermodynamic 
assumptions adopted as arbitrary conventions. The 
work of Guggenheim [7] as confirmed by de Groot 
and Tolhoek [8] shows that in ordinary electro- 
chemical systems only the term d^\ on the left in eq 
(1) may be included with cfln m\ as an operational 
quantity in a thermodynamic treatment. 

de Groot and Tolhoek give demonstrations that, 
in principle, one may treat transport across junctions 
between electrolytic solutions without seeking to 
make the division of chemical potential differentials 
called for in eq (1). The first object of this paper 
is to set down a fundamental steady state treatment 
which makes explicit use of familiar electrochemical 
quantities but avoids the use of eq (1). The funda- 
mental treatment given here is new only in the sense 
that it is a unified analysis applicable to the peculiar 
requirements of membrane studies. It may properly 
be viewed as first a restriction of Guggenheim's [9] 
treatment of galvanic cells to discontinuous systems. 
The treatment is then extended to nonisopiestic 
systems by including the flow of solvent in a mem- 
brane fixed frame of reference. It is further extended 
to the treatment of stationary states involving flows 
of electric charge. Finally, by following Temkin 
and Koroshin [10, 11] and Agar [12] in the treatment 
of electronic transport entropies we extend the 
analysis to nonisothermal systems. 

In much theoretical and experimental work there 
appears to be a need to adopt conventions regarding 
single ion ^ activities, membrane potentials, salt 
bridges, or ideal ''inert'' electrodes. In terms of the 
steady state theory these conventions are viewed as 
giving definitions to a ''force" — the EMF — conjugate 
to a "flow'' — the electric current — in the same sense 
that a pressure difference is a force conjugate to the 
flow of fluid volume. In the discussion given in sec- 
tion 10 we show an important case in which the ex- 
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plicit need for such a convention arises. The second 
object of this paper is to study the requirements to be 
met by an EMF convention if it is to lead to an in- 
ternally consistent treatment of transport phenom- 
ena. We will show that the application of the steady 
state theory leads quite naturally to a particular 
conventional definition of the membrane EMF. This 
definition is developed in sections 8 and 9. 

2. Analytical Model 

In our analysis we will formally treat a one dimen- 
sional, discontinuous, open steady state system with- 
out chemistry [13]. Flows through a junction are 
treated as if they took place between two points. 
The points, a and jS, are taken as representative of 
two parallel planes in a real system at each of which 
one may assign by experiments a truly representa- 
tive average value of any intensive property. A 
^'junction'* is everything between a and /? and its 
composition need not be defined. We only require 
flow continuity between a and jS. The flow of a con- 
served quantity, i.e., energy, mass as a chemical con- 
stituent, and electric charge, across a in the direction 
of /3 must be equal to the flow of that quantity across 
/3 toward the surroundings of jS. Flows are positive in 
the direction a to fi and are to be measured with 
reference to a and jS. 

The postulate of local equilibrium is applied at a 
and /3. A complete set of intensive equilibrium ther- 
modynamic properties can be defined by the local 
temperature, pressures, and component concentra- 
tions when a and ^ are at the same gravitational po- 
tential with no other external fields operating. Dur- 
ing a steady state these intensive properties are to be 
maintained constant by exchanges with the sur- 
roundings. 

We give a 'priori emphasis to certain consequences 
of this model. The uncertainty which arises in dis- 
cussions of transference numbers and diffusion coeffi- 
cients continues to receive attention in the literature 
[14, 15, 16]. An analysis of the model chosen gives 
permeability characteristics as integral characteris- 
tics of a junction in a cell fix frame of reference. Thus 
if the real system were a simple Hittorf transference 
cell [17], an integral transference number of an ion 
constituent would be by definition a property of the 
middle cell region, a solution filled glass tube. A 
transformation to an ordinary Hittorf or Washburn 
[12] number given with reference to the center of mass 
of the solvent would yield another integral property 
of the solution filled glass tube. The identification of 
these latter properties as properties of the solution 
would, in the usual way, require a demonstration 
that wall effects in large tubes contribute only to the 
extent that the tube defines the geometry of the 
solution. 

3. Fundamental Steady State Formulation 

de Groot gives a detailed discussion of the funda- 
mental principles of the steady state theory. We 
accept the theory here in the manner of a formal 



postulate. The material of this section sets forth the 
way in which quantities of established operational 
significance are to be employed to define flows and 
forces in the steady state description of junction 
processes without resort to eq (1). 

The choice of a concentration scale is immaterial 
to the essential results of this paper. In noniso- 
thermal systems it is generally preferable to avoid 
volume concentrations [37]. We find it convenient 
to employ molalities, mole kg~^ of solvent. The com- 
positions of the solutions at a and ^ are defined by 
the concentrations of r constituents. The first m, 
including the solvent, are nonelectrolytes. The next 
r-m are ionic constituents. We define a constituent 
composition vector, [m^], the elements of which are 
the molalities of these constituents. We define the 
constituent electro valence vector, [z^], where the 
elements are the electrovalences of the respective 
constituents represented in { m^} . The restriction of 
the electroneutrality condition is expressed by 



[z^Y[m']=0 



(3) 



where the superscript dagger, t> denotes transpo- 
sition. 

Under the restriction imposed by eq (3) it will 
always be possible to express the composition of solu- 
tions a and ^ in terms of r — 1 independent components 
[38]. We define a component composition vector, 
{m*}, where the respective elements are the molali- 
ties of electrically neutral isolable chemical com- 
pounds. These composition vectors must be con- 
nected by a relation, 



[^]{m*] = {m^ 



(4) 



or, in more detail, 
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An element, Vfjy of the matrix {v} is zero or a positive 
integer. It is the number of moles of the iili constit- 
uent contained in one mole of the isolated ^*th neu- 
tral component. The particular combinations of 
ionic constituents chosen as neutral components are 
a matter of convenience except that the independence 
of the components requires that {v} be of rank r—1. 
A unique solution must exist for the molalities of 
components in terms of molalities of r— 1 constituents. 
We note that the requirement of electroneutrality 
for every component is expressed by 



{znMHoi 



(6) 



We also note that convenience is best served if {v} is 
chosen such that the elements of the solution vector 
{m*j of eq (4j satisfy the condition 



m5>0; 



(i=l,2, 



-1). 



(7) 



It is important to note that in equilibrium thermo- 
dynamics one can treat the properties of a three 
component system such as H20-NaCl-KBr. In a 
steady state treatment of transport processes — 
'Vithout chemistry'^ — such as the interdiffusion of 
the components between solutions with differing 
concentrations of NaCl and KBr we must regard 
such solutions as special cases in whicli the concen- 
tration of a fourth component, either KCl or NaBr, 
is zero. In practice such special cases would be 
treated as the result of special restraints upon the 
possible transport processes in a system under study. 
Writing 



A/Xi=M?— M?. 



(8) 



as the difference in chemical potential of a chemical 
component or constituent between a and jS and 
using the equilibrium property 



M*=Z) ^ziM/i (i==l, . . ., r-1) 



(9) 
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(10) 



The superscript, *, indicates quantities which are 
functions of only temperatures, compositions, and 
pressures in accord with classical thermodynamics. 



The superscript, e, indicates quantities requiring an 
auxiliary electrochemical operation involving meas- 
ured transfers of electric charge coupled with 
measured transfers of material constituents [19, 20]. 
The algebraic theorem that eq (10) possesses no 
unique solution for the elements of the vector on the 
left is equivalent to Guggenlieiiu's conclusion that 
where a and ^ are formally treated as electrically 
isolated phases, the elements Ajuf have no general 
thermodynamic significance. However, when a and 
jS are not electrically isolated but comnmnicate 
through a junction permeable to at least one ion 
constituent, the elements A/x^ are operationally mean- 
ingful. Guggenheim considers the case when a and 
/3 are at the same temperature and pressure. He 
shows that when a reversible current can be passed 
between the terminals of a pair of identical electrodes 
reversible to the ith constituent, one at a and one at 
13, the relation 



AT.Pf^t=z{^(EI-E^) 



(11) 



defines the chemical potential differences when 
E^ — Ei is the electric potential difference measured 
between tlie electrode terminals. Equation (11) is 
easily modified to include cases where the tempera- 
tures and pressures at the electrodes are not equal. 
Here we follow Temkin and Koroshin [10, 11], 
Agar [12], and deBethune [21] by choosing to have 
the leads from the electrodes at a and /5 be wires of 
identical composition which come out to terminals 
which are at the same temperature, but we specif- 
ically require the use of platinum wire for these 
leads.2 We write tlie electrode reaction 



k ^i 



(12) 



for an electrode reversible to the ith constituent. 
jj^{ denotes the chemical symbol of the constituent 
corresponding to the ith element of {m^}- s/ik denotes 
neutral components of the electrode which are 
virtually insoluble in the solutions at the electrodes 
and, hence, are not represented in the vector {m*}. 
For example, for a silver-silver chloride electrode we 
write 



-Ag+AgCl-Cl- + 6-=0 



(12a) 



with k taking the values 1 and 2. s/n is silver; j3^2 is 
silver chloride; and X^i and Xj2 are —1 and 1, re- 
spectively. 

We then write for each ion constituent 



hlx\=z{^A<j>^] (i^m + l, . . ., r) 



(13) 



2 In nonisothermal systems it is necessary to indicate the reference medium 
implied for electrons when one writes an electrode reaction. The established use 
of platinum in constructing hydrogen electrodes and in resistance thermometers 
makes it convenient to regard platinum at 25 °C as the reference medium for 
electrons. 
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where A<?{>j is the EMF obtained by correcting the 
measured EMF, E^i—E", for the homogeneous 
thermoelectric effect in the platinum leads and for 
the chemical potential differences of the insoluble 
components involved in the electrode reaction, i.e., 

A.^,= (£■?-£?) +%Ar-^S X«Am. (14) 

where S'pt is the transported entropy for electrons in 

platinum.^ 

For silver-silver chloride electrodes under ordinary 

conditions eq (14) becomes 

-;|^[(Sl,-^A,ci)Ar-(FA.-FA,ci)AP] (14a) 

where the term in square brackets is derived from eq 
(16) below with the condition that 



(d<)^^''"-'^ 



(15 



for compounds in their standard states. 

Equation (13) defines the last r-m elements of the 
vector on the left in eq (10). The first m elements 
are defined using the classical equilibrium thermo- 
dynamic relation resulting from the choice of tem- 
perature, pressure, and molality as independent 
intensive variables, i.e., 

d^^,=-S'ldT^YUP^^{^!^ dm"^ (16) 

j=2 \Omj/T,P,m 

Si denotes the partial molal entropy; Vf denotes 
the partial molal volume; and the summation is over 
the independent solute components represented in 
the vector {m*}. Since we wish to express the 
thermodypamic ^^forces'', A/x*, as linear functions of 
the experimental ' ^forces' ^ it is convenient to write 



Afih 



—S^AT+VtAP+J2'%A^og m] 

3=2 



where 



e? 



VdlogmJ'/r.i 



(17) 



(18) 



with log denoting logarithms to the base 10. This 
gives a condensed symbol for our constant coefficient; 
we expect Ajit* to be more nearly proportional to 
Alog m* than to Am*. We restrict our attention 



8 At temperatures near 25 °C the work of Koroshin. and Temkin gives 
fpt=— 8.66X10-6+0.044X10-6T; (T>220 °k) 
-with temperature in degrees Kelvin and Spt in volt faraday deg-i. 



here to steady states representing displacements 
from a definite reference state of complete equi- 
librium between a and /S when the solutions a and jS 
have the same temperature, pressure, and con- 
centrations of components with no net flows of 
energy, chemical constituents, or electric charge 
taking place between them.^ The vectors {m^} and 
{m*} are defined at this state. 

It is useful to show how Oj is to be expressed in 
terms of the solution compositions and activity 
coefficients. 

We have the defining relation 



(Mf)r.P-M?+i?ri:^.anm{7{; 

k-\ 



(^-l,2,...,r-l)(19) 



where the summations are over the elements of {m^}. 
Here, for ion, constituents, 7^ is a single ion activity 
coefficient ^ which in this paper has no more than 
the ordinary significance of a formal device. It 
disappears when eq (19) is written in the meaningful 
form 



{iA)T,P=iA+RTiZvui In mi+vUTln 

k=l 



7? 



where 



v1='t^vu 



(20) 



(21) 



and 7* is the molal activity coefficient for nonelec- 
trolytes and the mean molal activity coefficient for 
electrolytes. From eq (20) and the definition 
given by eq (18) we find 



0* 



RT r , 
0.4343 L"^^ 



. Vk.VU : 



'mi 



-.r(l^) T (22) 



There is a real advantage of simplicity and con- 
venience in the form of eq (22) where ionic molalities 
are retained. The introduction of mean ionic 
molalities as suggested by textbooks leads to 
an almost hopelessly unwieldly expression in the 
general case. 

We recognize the experimental fact that the solu- 
tions may be of such a composition that no reversible 
electrode will function as required by eq (14). 
However, as demonstrated in section 6, if a single 
pair of reversible electrodes is available to measure 
A0i • for any one of the ion constituents represented 
in {m^}, the added equation supplements eq (10) to 
give a solution for the remaining terms, Afit Thus, 
given a satisfactory pair of reversible electrodes, we 
may define every term in eq (10) in terms of oper- 
ational quantities without resort to eq (1). 



i The steady state postulate can be applied with reference to other equilibrium 
states, e.g., osmotic equilibrium. JBowever, in writing eq (17) we imply the one 
indicated because we regard -Si*; Vi* and On* as constants. 

6 It should be noticed that we distingush 7* in eq (2) and 7/ here. We regard 
7« as being defined on the basis of a postulate that the electrostatic potential, 
A^, can be operationally defined between solutions of different compositions at 
finite concentrations. This postulate has been rejected. 
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with 



A fundamental set of phenomenological relations 

[j'}=im{xf} (23) 

IZl^IZ^t (24) 



may now be set down by choosing 



[J']- 



J{ 



JL 



Jm+i y , and {X^} = H 



Jr 



— AMm 

— AMm+1 

-Am^. 



(25) 



^J, J V.-Alnrj 

This choice of flows and forces is similar to that em- 
ployed by Kirkwood [22] in differential foiin. A\e 
have already defined the first r elements of {X^} 
and grant that AlnT is a meaningful quantity. The 
first r elements of {J^} are flows of chemical con- 
stituents corresponding to the respective elements of 
{m^}. Since the chemical constituents are defined 
as having fixed identities in cell processes including 
exchanges with the surroundings of a and /? it may 
be granted that the first r elements of {J^} are 
meaningful quantities. It is this conservation of 
chemical identity in the complete specification of 
concentrations and of flows of matter which permits 
the analysis ' Vithout chemistry" [13].^ 

The flow of '^entropic heat," JJ, must be given 
special notice. We select the name, entropic heat, 
on the basis of the discussion set down by de Groot 
[24] who notes that various flows are termed ^^flows 
of heat" by different authors. The flow of entropic 
heat is defined by 



Ji=ju-'i:i^ij{- 2 Mfj{ (26) 

1=1 i=m-f-l 

where J„ is the flow of energy. Energy is of course 
conserved in all cell processes. Entropic heat is not. 
However, we cannot rely upon eq (26) alone to 
define J^ since the chemical potential of a single ion 
constituent, ^f, has been given no meaningful defini- 
tion. Only its significance in the particular linear 
combinations of eq (9) has been established. It is 
necessary to state conditions which are at least 
sufficient to establish an operational definition of JJ. 



We first note that by eq (23) 

JJ- -L J,A InT- 1] LijX^ 



;=i 



(27) 



6 The treatment of water as a nonelectrolytic solvent limits th'^ participation 
of hydrogen and hydroxyl ions in the net transfer of electric charge to those sys- 
tems where either a strong acid or strong base is present as a solute. The treat- 
ment of weak electrolytes requires the introduction of "chemistry" [23]. 



A measurement of the thermal conductivity of the 
junction when all flows of matter are zero gives 



Ji=Ju', (Jg.=0) (28) 

from which i^^ may be determined using 



LL 



^" '-±LfjT^'^ (Jj^a=0) (29) 



AlnT ' pi '^ AlnT 





U,=u, 


equations 




Li,=U,= 


AlnT' ^i "AlnT 



when the experiment is designed to simultaneously 
evaluate the elements of {X^} and the elements 
L^j are already known. Since by eq (24) 

(30) 



(31) 



applied to measurements of the flow of each material 
constituent in the presence of a temperature differ- 
ence would evaluate each term L^gj when again the 
elements of {X^} are simultaneously evaluated and 
the elements Lij are known. We finally require the 
measurement of a complete set of isothermal per- 
meability characteristics to define the elements of 
the isothermal admittance matrix in the equation 



■LU 



J f 



iLL'rl ■ ■ ■ Ur. 



UYd \ju 



(32) 



which is obtained from eq (23) by setting AlnT 
equal zero. Isothermal measurements dealing with 
flows of conserved quantities require no special 
discussion. 

This formulation of a linear macroscopic theory 
requires an additional postulate to put a restriction 
upon the magnitude of any admissible displacement 
of equilibrium from the reference state. We have 
chosen to give the treatment of a discontinuous 
system in order to avoid operational uncertainties 
which often arise in attempts to apply an analysis 
of a continuous system [25]. In general we are 
unable to measure intensive properties at an arbi- 
trarily dense succession of planes lying between 
a and fi. Our choice is in line with the suggestions 
of Kirkwood [22] and Scatchard [26]. Thus the 
specification of a junction such as one including a 
membrane or porous plug must include the specifi- 
cation of the composition of the contacting solutions 
at the reference state. At the reference state we 
formally assume that equilibrium is established with 
respect to all solution components throughout the 
region a to /?. The integral admittance matrix, 
\L% can in principle include linear ^'polarization'' 
effects in the boundary layers. However, with 
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reference to Kirkwood's treatment of thed tirec 
transformation of the analysis of a continuous 
system into the analysis of a discontinuous system, 
Schlogl [27] has demonstrated the possibility of 
^'apparently'' linear phenomenological relations 
where the elements of {Uj are not proper constants. 
In such a case the respective elements of a series of 
differential admittance matrices defined as functions 
of position between a and /3 would not, upon inte- 
gration, give the corresponding elements of [il. 
In a purely phenomenological approach no basis 
can exist for distinguishing ' 'really'' and ''appar- 
ently" valid applications of eqs (23) and (24). 
Schlogl's analysis suggests that no difficulty exists 
in the limit of the reference state and that if the 
forces represented in {X^} are "sufficiently small," 
the experimental demonstration of any such dis- 
tinction would be difficult and would require meas- 
urements of the highest accuracy and precision. 
We require some quantitative criterion for judging 
what forces are "sufficiently small." 

Miller [28] has recently reviewed the experimental 
evidence which bears upon the validity of the 
Onsager relations in macroscopic processes. Miller 
[29] and Dunlop and Gosting [30, 31] have givea 
detailed attention to ternary diffusion in electro- 
lytic solutions. It is clear from a study of these 
papers that any criterion will necessarily be some- 
what arbitrary. Nevertheless, since a major part 
of the intuitive support for the postulation of a 
linear macroscopic theory involves the established 
validity of linear free energy relationships such as 
eq (17) it is reasonable to assume that the range of 
validity of eq (17) is of fundamental significance in 
applications of the theory. Therefore, we choose to 
apply the following minimum condition : 

It is assumed that determinations of the elements 
of {Uj will involve displacements of equilibrium 
with forces and flows of such a magnitude that the 
uncertainty in the value of an element L{j arising 
from uncertainties in measurements of temperatures, 
concentrations, pressures, and potentials will be 
equal to or greater than the errors introduced by 
neglecting variations of F*,>S*, and 6*^. 

This criterion has the practical advantage of 
setting "sufl&ciently small" at magnitudes consistent 
with the precision and accuracy of any particular 
experimental investigation. 



4. Transformation Properties 

In the fundamental formulation, a flow of electric 
charge is identified with the flow of each ionic con- 
stituent and a conjugate electromotive force is de- 
fined for each such flow. However, no net flow of 
charge or overall electromotive force occurs as an 
element of the vectors [J^] and {X^}. The thermo- 
dynamic theory does not require us to define vectors 
containing such elements, but the definition of all 



elements of {J^} during all processes including 
electric charge flows depends upon our extra- 
thermodynamic definitions of ionic constituents 
and our use of the established conservation relation 



^ i= m-1-1 



(33) 



The electric current in amperes is denoted by /. 
Equation (33) suggests a linear transformation of 
flows to give a new flow vector, say {«/'}, where one 
element is the flow of electric charge. A correspond- 
ing transformation of forces to give a new force 
vector, {X'}, would contain an element conjugate 
to the flow of charge. In this paper we identify 
such an element of {X'} as a junction EMF. 

The steady state theory explicitly deflnes restric- 
tive conditions which must be applied in carrying 
out linear transformations of flows and forces. The 
transformation properties are summarized by de 
Groot [32]. We set them down briefly in the nota- 
tion of this paper. 

We take as given 



with 
and 



lLq=iuf 

Ta={JfY{Xr} 



(23) 
(24) 
(34) 



where cr is the rate of entropy production during a 
steady state process described by the vectors [J^] 
and {X^}. 

We may choose an alternative description of the 
same steady state process in terms of vectors { J^} 
and [X^] where 





{J']=IL'Y[X'} 


(35) 


with 








IL'l^Wf 


(36) 


and 








r<r={J'}t{Z'} 


(37 



The entropy production, cr, must be invariant under 
a linear transformation of flows and forces. We 
must define a new force vector by 



{X'}=W'l{Xf} 



(38) 



where the order of the superscripts in the matrix of 
the linear transformation indicates the direction 
of the transformation, i.e., 



wi=W'r' 



(39) 



The elements of the vector {X'] must be hnearly 
independent combinations of the elements of {X'], 

detWly^Q (40) 
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where the prefix, det, denotes the determinant of 
tlie matrix. The corresponding transformation of 
flows is uniquely defined by 



where 






(41) 
(42) 



Alternatively, since the reciprocal of the transpose 
of a matrix is the same as the transpose of its reciprocal 
the transformations may be defined in reverse order 



Mntli 



/'^ll = fra>'']lt-l 



(43) 



The matrix fa^'] takes the role of de Groot's matrix 
jS and we note that an application of the trans- 
tormation 



lUl=W^\lUlia- 



,^'ilt 



(44) 



is required to complete the definition of an alternative 
set of phenomenological relations to replace eq (23). 
Finally, we note that in the general case we are 
dealing with a force vector whose elements are 
formally linearl^^ independent in terms of their 
definitions from measurements of intensive prop- 
erties. We may state that if not all the elemeuls 
X{ are zero. 



2 6,X;{^o 



;=i, 



(45) 



for any choice of constants, ei. 'J'he trausformation 
properties given remain valid regardless of the 
permeability characteristics of any junction perme- 
able to at least one ion constituent. In a real 
system if any one flow or a linear combination of 
flows is identically zero in all steady state flow 
processes, the elements of {J^} will not be linearly 
independent. We will study such a special case 
in section 7 below and set a linear combination of 
forces to zero, de Groot demonstrates that the 
a priori definition of independent forces secures the 
validity of the transformation properties despite 
any linear dependence of flows in a special case. 

5. Restrictive Condition Upon a Defined 
EMF 

We consider first a matrix [a^'I- Certain general 
restrictions upon the form, of this matrix may be 
adopted on the basis of elementar^^ considerations. 
The non electrolytes are independent com.ponents of 
conserved chemical identity. The law of conserva- 
tion of m.ass is applied to the flow of each such 
nonelectrolyte. The flows of r—m—1 electrolytes 
and of electric charge are an interdependent group 
subject to the laws of conservation of mass and of 
electric charge, but the electroneutrality condition 
only operates as a restriction upon the flow of ionic 
constituents with respect to other ionic constituents. 
It does not restrict the independence between flows 



of nonelectrolytes and electrol.ytes. We take these 
points into consideration when we write 



1 . 


. 


1 




. 










. 


. . 1 


1 





. 


. 


0-m+l 


m-fl • • 


d'm + l, T- 


1 (lm+\,r 





. 


. 


'«r-l, 


w+1 • • 


(lr-i,T-l 


dr-hr 





. 


. 


' 7f 

1 ''m+l 






z{ 







1 




. 







. 


. 


1 



(46) 



where the last column and the last row indicate that 
we retain the flow of entropic heat as a fundamental 
flow. It is convenient to have the rth element of 
our new vector {J'} equal to the flow of electric cur- 
rent. This is indicated in tlie matrix where 



and 



a{;=0; (i-1, . . ., m,g); (47) 

<=2.^; (i=m+l,...,r). (48) 



The rth flow defines the flow of electric charge. 
The flows J[ through J'm are the flows of neutral 
components not involving a transfer of electric 
charge. It is natural to inquire into the possibihty 
of defining as set of r — m—1 independent flows of the 
neutral electrolytic components represented by tlie 
last r—m—1 elements of the vector {m*} in such a 
way that these flows do not imply a flow of electric 
cliargc. On the most elementary rational grounds 
a net transfer of a quantity of an electrolytic com- 
ponent from aio ^ must represent some combination 
of the transfers of ionic constituents. It is not 
obvious that any established rules or conventions 
dictate these combinations, but we identify the 
respective conjugate forces defined by a matrix 
f6-^'J with the chemical potential differences of 
neutral electrotytic components defined as linear 
combinations in eq (10). The fact that the eq (42) 
describes a unique operation permits us to carry 
out an inverted development. 

Before setting down a matrix {¥'} we note that in 
tlie m.atrix fa^'J given above the dashed lines show 
how the matrix may be partitioned to give a pseudo- 
diagonal matrix; 



K'I= 
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^J 



(49) 



The submatrices filf f and IQ\ are identity matrices, 
and the submatrices {01 are zero matrices. It re- 
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suits as a property of such forms [33] that when the 
matrix {N} is not singular, 



lbf'}=wv-'= 



with 



M 0" 

A^t-iQ 
LO QJ 



Thus we write the pseudodiagonal form 



(50) 
(51) 
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(52) 

with the matrix \y\\ from eq (10) comprising the 
first r— 1 rows and first r columns. Only the ele- 
ments, h{'j, in the rth row and columns m+1 through 
r are left to be determined. 

The restriction upon IV} can be demonstrated. 
For the moment we assume that det \¥'\ does not 
vanish. The definition of an element of the ith row 
of the reciprocal of a transposed matrix is [33], 






(53) 



where the numerator on the right is the cofactor of 
the element indicated. Following eq (53) we write 
the definitions of the elements of W\ within the 
submatrix {N} as 



ff COJ -^^ i-m,j-rri 



(54) 



The electrical neutrality of each component repre- 



sented in the vector 
the condition, 

^ I ^j ■'■^ i-m,j-m ^} 
j=m-{-l 



m*} as stated in eq (6) gives 



(^=m + l,...,r-l) (55) 



Equation (55) applied in the first r—m—1 rows of 
det IN}^~'^ gives zeros as the first r—m—1 elements 
of any column when the other r—m—1 columns are 
added to it. The last element of that column 
then becomes 



S ^f iv;-v,--.- S 2^H'j (56) 

j = m-\-l j=w+l 



for any jth column. It then follows that 
detW'}='-^( ± z^y/) 



(57) 



for any jth column m-fl through r. Substituting 
from eq (57) for det {¥'} in eq (53) gives 






Clearly, if 



^/(i: ^/&^;) ' (58) 

\j=m-hl / 



j=m-i-l 



(59) 



eqs (47) and (48) will be satisfied provided (¥'} is 
not singular. The nonsingularity of {b^'j is assured 
by the requirement that {v} be of rank r— 1. At 
least one cof bi] will not vanish. When the condition 
of eq (59) is applied the determinant will be given by 



6. Reference Ion Electromotive Force 



(60) 



An extremely important class of experimental 
systems including concentration cells and various 
cells with transference have identical reversible 
terminal electrodes [1] and can be designed to fall 
into the class of systems discussed ia section 2. It 
is proper to inquire as to whether or not the EMF 
between a single pair of identical reversible probe 
electrodes, one at 4> and one at fiy may be formally 
regarded as the junction electromotive force. We 
may choose to order the elements of {m^}, {X^} and 
{J^} to have the rth constituent be that one to which 
the electrode pair is reversible. We define a ''refer- 
ence ion" force vector with 



Xl=X{] (^-l,...,m); 



(61) 



X[= S ^.z^/; (i-m + l,...,r-l); (62) 



^' z' 



and 



XI 



-X{- 



(63) 



(64) 



The matrix [¥''] is thus of the class of [¥'] in eq (52) . 
The elements 



and 



H'i=^; {j=l,...,r-l,q) 



^ TT 



(65) 
(66) 



clearly satisfy the restriction imposed by eq (59). 
Therefore, any one of the r—m EMF's A^^ defined by 
eq (14) may be chosen as the junction EMF. 

This reference ion formulation is important for 
practical applications. In most systems one can at 
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best discover a single pair of adequate reversible 
electrodes. Experimental flows will normally be 
expressed most directly in terms of the vector { J^} 
with the forces being in terms of the vector {X^}, 
A transformation, 



must be applied with 



(67) 



(68) 



7. Junctions With Zero Electric Current Flow 

Treatments of electromotive forces acting across 
junctions are most often encountered in discussions 
of the thermodynamics of galvanic cells. In such 
treatments one is concerned with the apparently re- 
versible EMF measured at the lunit of zero electric 
current flow. Guggenheim [9] has given a treatment 
of isothermal, isopiestic cells with liquid-liquid junc- 
tions without employing eq (1). The treatment is 
one of a continuous steady state system without 
chemistry in a solvent fixed frame of reference. His 
work is an exception to the apparent universal prac- 
tice of using eq (1) or its equivalent. Therefore, 
some indication must be given of the relationship 
between the treatment here and the more common 
treatments. In addition the treatment of junctions 
at states of zero current gives an important relation- 
ship for use in connection with eq (10). 

It can be shown [34] that when no electric current 
is flowing through the junction, 



Sr/A/x.-O; 



(i=»). 



(69) 



The coefficients, t{, are the integral stoichiometric 
transference numbers in a cell fixed frame of reference 
when no difference of temperature, pressure, or com- 
position exists between a and jS, i.e., 

ri^jdj: 2'j jiy\ (Ar-AP-{Am*}=0). 

(70) 

The numbers r/ are properties of the junction defined 
as functions of the elements of the admittance ma- 
trix iiq. 

These functions are 



1 



where 



with 



■Lie 1 = 1 



j=l i=l 






(71) 



(72) 



(73) 



The detailed manipulations leading to the above 
relations need not be repeated. However, in existing 
analyses the use of eq (1) as an elementary postulate 
gives a familiar result in differential form as 



-^d^=^r^j(dfJijU=o=RTj:r^j{dln a,),^=o (74) 

7=1 j=l 

which because of uncertainties in notation may ap- 
pear to be inconsistent with eq (69). The relation- 
ship may be seen clearly if we write 



A/Xj=(AhjU^o+^'j^^^ 



(75) 



and substitute for Ajuy in eq (69). The use of eq (73) 
to reduce the coefficient of J^^ then gives 



-J^A^ 






(76) 



Of course this only demonstrates a formal abstract 
relationship since the operational significance of eq 
(1) has not been established for ordinary electro- 
chemical systems. 

Nevertheless, we note that eq (69) is a hnear com- 
bination of chemical potential differences defined as a 
function of the permeability characteristics of the 
junction. Thus when this relation is added to the set 
of equations in eq (10) we have a set of r equations in 
r unknowns, 
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y (77) 



The elements r^+i through t{. obey eq (73) ; there- 
fore the condition of eq (59) required for J^'^+i 
through ¥/r ^^ treating the matrix {¥'} applies. It 
may be shown by a treatment similar to that applied 
to {¥'} that the matrix here does not vanish unless 



n+l- 



=0 



(78) 
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Thus if the junction is permeable to at least one ionic 
constituent of the solutions at a and /3, eq (77) pos- 
sesses a unique solution. This result confirms the 
statement given after eq (10). 

It should be noted that eq (77) has the distinct 
practical advantage of unifying the treatment of 
galvanic cells having identical terminal electrodes in 
compartments separated by a defined junction 
through which flow continuity is established. It is 
clear that a solution for any A/xf will after applying 
eq (13) give the EMF, A0,, at the limit of zero electric 
current flow in terms of the transference numbers 
and the thermodynamic properties of independent 
isolable components. 

8. Electromotive Force at Infinite Dilution 

The reference state of unit activity coefficient es- 
tablished for treating the thermodynamic properties 
of solutes is infinite dilution in pure solvent. In the 
limit of infinite dilution aqueous solutions of strong 
electrolytes are ^^nearly" insulators and the Debye- 
Htickel limiting law is valid for electrolytic solutes as 
a whole. In accord with the discussion of Tolhoek 
and de Groot [81 and Guggenheim [7] we may regard 
it as meaningful to speak of a difference of electro- 
static potential acting as a part of the ^^force'' con- 
jugate to the transfer of an ionic constituent between 
solutions of infinitesimally different compositions in 
the neighborhood of infinite dilution. To define this 
'^force'' we make use of several commonly used extra- 
thermodynamic conventions, and we need only con- 
sider isothermal, isopiestic systems. 

We assume that the activity coefficients of solute 
nonelectrolytes approach unity much more rapidly 
than their concentrations approach zero in the limit 
as a solution is diluted with pure solvent. We write 



(rfM*)2 



=^rainm?; (?;-2, . . .,m). (79) 



For ionic solute constituents we apply the same 
assumption and write 

W^T,p=^nTd\^m{^z{3^di) (i-m+1,. ..,r) (80) 

It will be noted that the Debye-Hiickel theory sup- 
ports the assumption underlying eq (80) since the 
limiting law 

(ln7l)T.P..^=o— ^(4)'V^ (81) 



where u is the ionic strength, i.e., 



(82) 



suggests that (^\iY^T.p,di=^ approaches zero in the 
limit of infinite dilution as In m{ approaches minus 
infinity. 

Using the relation. 

d\ixm{=^ (83) 

in eqs (79) and (80) and then applying the electro- 
neutralitv condition 



2 2{dm{=0 

;=m+l 

gives a set of r equations in r unknowns. 



(84) 
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Tiie coefficient matrix is nonsingular; its determi- 
nant is given by 



2{BT)'-'M(n m{\ >0. 



The solution, 



S^d^-- 



j=m-\-l 



4? ('^''^)-.-- 



(86) 



(87) 



gives the differential of the electrostatic potential at 
infinite dilution as a linear combination of the 
differentials of chemical potential. The coeffi- 
cients in eq (87) satisfy eq (59). 

It is clear from the form of eq (80) that d^ is 
introduced as the '^force^' conjugate to any infini- 
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tesimal transfer of electric charge, hence, eq (87) 
defines the electromotive force at infinite dilution. 
We use eq (87) to define the elements of the rth 
row of the linear transformation matrix in 



l¥q{Xn = {X'} 



(88) 



where l¥^ is of the class required by eq (52). By 
eqs (41) and (43) we have the transformation 



lafq{Jf} = {J<^} 



(89) 



where each element J^+i through J^-i is a linear 
combination of the flows of ionic constituents de- 
fining the '^flow'' conjugate to the ^'force" acting 
on an isolable electrolytic component as a whole. 
The combinations are necessary consequences of 
established extra-thermodynamic conventions and 
assumptions. They are implicit in our established 
methods of treating the properties of completely 
dissociated electrolytes. 

9. The Membrane EMF at Higher Concen- 
trations 

The definition of the membrane EMF at higher 
concentrations follows directly from the definition 
at infinite dilution. In order to remove any implica- 
tion of a definition of an electrostatic potential 
difference at finite concentrations and concentration 
differences we write 



z{m{ 



i= OT+1 -iU 



Am? 



(90) 



where ^ is the electromotive force between a and /?. 
We regard \[/ as representing the action of an electric 
field in the sense that it denotes a force acting to 
produce a flow of electric charge, but it is not a 
!:> gradient of an external electrostatic potential field 
i [35]. The EMF, xj/, is fully meaningful since every 
|) term on the right hand side is meaningful. It is 
I arbitrary in the sense that our choice of the reference 
r state of infinite dilution and the applications of eqs 
(79) and (80) are arbitrary. It depends upon our 
definitions of the elements of the vector {m^} and 
^ our use of these ion constituent concentrations in 
writing the ionic strength, u, given in eq (82) and 
p the electroneutrality condition given as eq (3). 
! It is, however, quite independent of any arbitrary 
[' choice of the matrix {v} to give the elements of 
I {m*}. The matrix {a^^ of eq (89) defined at infinite 
I dilution fixes the stoichiometry governing exchanges 
I of radical chemical constituents. The most elemen- 
;^ tary considerations of internal consistency dictate 
I the use of this same stoichiometry regardless of the 
[ absolute concentrations of the solutions. 

I 10. Discussion 

j^ The need for defining a membrane electromotive 
{ force and the corollary flows of electrolytes as n. 
y whole arises quite naturally if we consider the treat- 



ment of electrokinetic phenomena. However, we 
must first establish certain points in terms of a 
specific example. 

Consider a membrane cell with conipartnients 
filled with HCl solutions. Let it be arranged in 
the manner of a four lead resistor witli a pair of 
probe electrodes close to the membrane — one on 
each side. These probe electrodes establisli the 
positions of the planes, a and P, indicated in section 2. 
They are potential indicating devices and carry no 
steady electric currents. Let each compartment 
also be provided with a working or current electrode 
more remote from the membrane than a and /3. 
The electrode reactions at the current electrodes 
are to involve virtually insoluble components and 
either the chloride or the hydrogen ion constituent 
of the solution. Let us also provide for continuous, 
adjustable flows through each compartment utilizing 
feed solutions of adjustable composition. 

A little reflection will show that we can establish 
a particular steady state — defined by giving the 
elements of {X^} — regardless of whether it is the 
hydrogen ions or the chloride ions which are ex- 
changed with the solutions at the current electrodes. 
The adjustments of feed solution compositions and 
rates of flow serve as our compensating variables. 
If the phenomenological relations, eq (23), describe 
a steady state process taking place between a and P, 
they describe it independently of the specific nature 
of the processes at current electrodes which may in 
practice be at virtually infinite distances from the 
region between a and /3. This statement of inde- 
pendence, the independence of dissipativc proc- 
esses occuring in different regions of space, expresses 
an elementary requirement of the steady state 
theory [39]. Processes which are described as 
coupled — liaving phenomenological relations with 
cross coefficients — in the production of entropy 
must take place in the same system. By system 
we mean an open system, the region of space be- 
tween a and /3 as defined in section 2. If we wish 
to describe a particular process in terms of some set 
of phenomenological relations resulting from a 
linear transformation of eq (23), the second form of 
the phenomenological relations will describe the 
same and only the same process if and only if it is 
defined so as not to include external dissipative 
processes. This elementary requirement is ex- 
pressed in the steady state theory when we state 
that the rate of entropy production given by eq (34) 
must be independent of a linear transformation of 
the forces and flows. 

Ordinary isothermal electrokinetic phenomena at 
membranes are investigated by making experimental 
measurements — in a membrane fixed frame of ref- 
erence — of flows of fluid volume and electric current 
under the influence of their conjugate forces — pres- 
sure difference and electric potential difference. 
We can site, for example, de Groot's [40] and Gug- 
genheim's [36] demonstrations of the application of 
the steady state theory in the special case where no 
differences of temperature or of composition are 
present to act as forces. They show that the well 
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known experimental relation, the Saxen relation, be- 
tween the streaming potential and the electroosmotic 
transport, i.e., 






m 



AP=Ar={Am* 



(91) 



may be regarded as a consequence of the Onsager 
reciprocal relations. JV is the rate of flow of fluid 
volume and A^ is, as before, the difference in elec- 
trostatic potential. The electrostatic potential, Af, 
is of course a meaningful quantity when measured 
between phases of identical composition and tem- 
perature; it is the EMF measured between the ter- 
minals of any pair of identical electrodes. Hence, 
neither of the treatments cited makes an}^ reference 
to the nature of the electrodes employed to measure 
AJ. The fact that these analyses impose the restric- 
tion 

Ar={Am*}=0 (92) 

and the fact that in practice the experimental 
measurements do not include measurements of the 
flows of heat and of the relative flows of constituents 
obscures certain features of the phenomena involved. 
Jv is computed from the measured rate of change of 
the volume of solution in a cell compartment with 
due allowance for the change of the volume of a 
compartment arising from the reactions of insoluble 
components at the working electrodes. In such 
treatments there is normally an implicit under- 
standing that the reactions at the probe and work- 
ing electrodes are identical. The problem of sep- 
arating electrode and membrane processes simply 
does not arise in such a restricted case. However, 
we may, for example, study Guggenheim's [9] treat- 
ment of concentration cells with transference. This 
is a case where {Am*} is not restricted to zero. Al- 
though our formulation is first to be distinguished 
from his treatment of concentration cells by the fact 
that his is in a Hittorf, i.e., solvent fixed, frame of 
reference, the essential point is that Guggenheim 
explicitly requires the use of electrodes reversible 
to one of the ionic constituents. Such a treatment 
avoids the classical uncertainty involving A^. 
Since he restricted his treatment to steady states 
not involving flows of electric charge, the require- 
ment of specific electrodes does not introduce any 
difficulty with respect to the separation of elec- 
trode and junction processes. 

Since the phenomenological relations of eq (23) 
describe processes which involve electrokinetic effects 
as well as diffusion effects, we are led immediately to 
attempt to set down a treatment in which heat flows 
and temperature differences and diffusion flows and 
concentration differences are added to the flows and 
forces of the electrokinetic treatment. However, 
we find that not only the meaning of AJ but also the 
meaning of JV, is uncertain. We must require our 
formulation to describe the processes in the ^ ^mem- 
brane'', i.e., a to /3, region. It is to result from a 
linear transformation of the fundamental phenom- 
enological relations, eq (23). 



At first glance it appears that we would write 

m r 

J'y=i:V(Jf+ S ViJ{ (93) 

where the terms V{ are the partial ^olal volumes of 
the nonelectrolytes and the terms F? are the partial 
molal volumes of ionic constituents represented in 
the vector {m^}. However, we have already indi- 
cated in the discussion of eq (2) that Ff has no opera- 
tional meaning when we reject the postulate that 
A^ is operationally defined. We could consider a 
convention which defines Ff for any constituent. 
Such a convention would permit us to retain the 
distinction between electrode and membrane proc- 
esses. However, there is an alternative approach 
which has to be considered. 

We recognize that in the study of combined 
hydrodynamic and diffusion flow processes diu*ing 
steady states not involving net flows of electric 
current, we have by definition 



r-l 



(94) 



where («/*)/=o are the flows of the neutralcomponents 
represented in the vector {m*} and Vf are their 
partial molal volumes. We recognize that in prac- 
tice the flows (J*)7=o will be computed from experi- 
mental measurements of changes of ion constituent 
concentrations with explicit dependence upon the 
requirement of electron eutrality, i.e., as linear com- 
binations of the flows (Ji)i=o. We are thus led to 
write 

J'y=^VJJ'j',l7^0 (95) 

where the flows Jj are linear combinations of the 
flows of ion constituents which represent flows of 
electrolytes as a whole. It was necessary to dis- 
cover what conditions are to be met in defining Jj. 

In section 5 of this paper we showed that we are 
first led to the conclusion that when we define the 
membrane EMF in terms of the potentials measured 
between pairs of identical probe electrodes with 
specific electrode reactions, we at the same time 
define, by implication, a particular method of com- 
bining the flows of ionic constituents into the flows 
of neutral electrolytes as a whole. Thus, by im- 
plication, we establish the definitions of Jj and, 
hence, Jv in eq (95). Therefore, we conclude that 
it is not proper _to introduce any auxiliary convention 
which defines Ff in eq (93) when we set down a 
transformation of eq (23) into an electrokinetic 
formulation. 

In section 6 we showed that we can choose a 
^ ^reference ion" EMF, A<^„ by choosing to employ 
probe electrodes reversible to some particular ion 
constituent. Formally, there will be r-m distinct 
choices possible. In general there will be a family 
of r-m pairs of matrices {¥''} and [a^''}. The defini- 
tion of the flow of any ^'th neutral electrolyte 
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j;= X: <J{;(i=m + l,...,r-l) (96) 

will in general be different for each choice of a 
* 'reference ion.'' Although the arbitrariness of a 
'^reference ion'' electromotive force is objectionable, 
one can always employ such a device in order to 
define electrokinetic effects in terms of a general 
treatment. We can employ the essentially ''direct 
experimental variables" such as appear in Saxen's 
relation, eq (91), in conjunction with diffusion flows 
and concentration difference forces without including 
external dissipative processes in our description of 
membrane processes. 

The arbitrariness of a "reference ion" treatment 
can only be removed by adopting an additional 
convention- — an "averaging" convention. In Sec- 
tions 8 and 9 we showed that such a convention is 
implicit in the conventions already established for 
the treatment of the thermodynamic properties of 
electrolytic solutions. A general treatment which 
is consistent with the elementary definmg stoichio- 
metric relationships employed at the reference state 
of infinite dilution results when we write 






(97) 



where Jj is defined by the transformation given in 
eq (89). 

It is important to emphasize that the treatment 
which defines \l/ is complete only with respect to the 
requirement of consistency with established conven- 
tion. We have implied that i/' is a physically sig- 
nificant quantity. The developments of this paper 
are not adequate to demonstrate any real physical 
significance. Our development in sections 8 and 9 
follows the formal procedure of defining the Lagran- 
gian multiplier corresponding to the electroneutral- 
ity condition.^ Although as a formal device in sec- 
tion 8x1/ is uniquely defined when we write eq (85), 
the approach to ideal behavior implied by eq (80) 
is not subject to an adequate experimental test [42]. 
The use of the concept of "almost" insulators is not 
satisfying as a basis for attaching physical signifi- 
cance to A^ in extremely dilute solution and, hence, 
to xp in more concentrated solutions. 

An important problem for further investigation is 
the comparison of experimentally determined values 
of \l/ with the EMF's measured between pairs of 
saturated calomel electrodes. A study of represen- 
tative attempts to set down meaningful conventions 
for junctions potentials [4, 43, 44] based upon estab- 
lished common practice shows that saturated solu- 
tions of potassium chloride are regarded as special 
experimental devices. It is abundantly clear that 
they are assumed to represent a practical approach 
to some ideal of "inert" behavior in experimental 
systems. We can, for example, accept t/^ as a phys- 
ically significant quantity if a general study of a 
large class of experimental systems shows that sat- 

7 Lagrangian multipliers as physically significant quantities are common in 
statistical mechanics. Onsager [41] suggests the use of such a multiplier to 
represent the electric potential in electrochemical systems. However, he gives 
no explicit definition and, further, he chooses to have his dissipation function 
include electrode processes. 



urated calomel electrodes can be regarded as exper- 
imental devices which measure xp directly with a 
degree of accuracy which is adequate for most prac- 
tical purposes — ^^to within a few millivolts." We 
suggest this hypothesis here because our first tests 
support it, but we regard it, for the present, as 
only an example. 

11. Summary of Symbols and Notation 
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Activity of a constituent, eq (74). 

Matrix of transformation of g= 
r-j-1 flows, r of chemical con- 
stituents and one of entropic 
heat, into the flows of r— 1 
chemical components, electric 
current, and entropic heat, eq 
(46). 

la^'Y], the reciprocal of fa^'J. 

Matrix of the class \a^'\ consistent 
with the conventional EMF, 
^, section 9. 

Matrix of the class fa-'''} consistent 
with a particular reference ion 
EMF, A</>„ section 6. 

Chemical symbol of electrode com- 
ponent, eq (12). 

Matrix of transformation of forces, 
eqs (50) and {b2)]l¥'}=^Wf-\ 

Matrix of the class {V'}] {¥^1= 

Matrix of the class W'h {¥'}= 

WT\ 

Chemical symbol of ion consti- 
tuent, eq (12). 

Cofactor of i/th element of \a/'\ 

Determinant of matrix fa^'J. 

EMF measured between electrode 
terminals. 

Faraday. 

Electric current. 

'1=1, . . ., r; rate of flow of a 
chemical constituent. 

%=l, . . ., r— 1; rate of flow of a 
neutral component. 

//J^eq (33). 

^=1, . . ., r— 1; rate of flow of a 
neutral component when 7=0, 
eq (94). 

g=r+l; rate of flow of entropic 
heat. 

Rate of flow of fluid volume, 
section 10. 

Admittance matrices. 

Constituent composition vector 
having r elements which are the 
molalities of m nonelectrolytes 
and r—m ion constituents. 

Component composition vector 
having r— 1 elements which are 
the molalities of m nonelec- 
trolytes and r—m—1 neutral 
electrolytes. 

Pressure. 
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B 

S*t 

-Si 

T_ 

Vf 

y\ 

XijX'i^XijXl 

X{ 

XijXijXi 

xi 

X'r 

X^qjX'q,Xg,X'^ 



2{ 
a,0 



7? 



7! 



7f 



Am? 



M? 



Gas constant. 

Partial molal entropy of iih com- 
ponent in {m*}. 

(Undefined) partial molal entropy 
of an ion constituent. 

Temperature. 

Partial molal volume of 7th com- 
ponent in {m*}. 

(Undefined) partial molal volume 
of an ion constituent. 

i=l, . . .J m; forces, —Am*, acting 
conjugate to the flows of non- 
electrolytes. 

i=m-{-\, . . ., r; forces, — Ajuf, 
acting conjugate to the flows of 
ion constituents. 

i=m+\, . . ., r—1; forces, —Am*, 
acting conjugate to the flows of 
neutral electrolytes. 

Conventionally averaged force, 
— J^i/', acting conjugate to the 
flow of electric charge, section 9. 

Force, — .#A<^;., acting conjugate to 
the flow of electric charge in 
accord with a reference ion 
convention, section 6. 

q=r-\~l] force, —A In T, acting 
conjugate to the flow of entropic 
heat. 

Electro valence of a constituent. 

Denote boundaries of membrane 
region; a to P positive. 

i= 1 , . . . , m ; molal activity coeflS- 
cient of a nonelectrolyte; and 
i=m-yi . . ., r— 1; mean molal 
activity coefficient of an elec- 
trolyte. 

(Undefined) electrostatic molal 
activity coeflScient of an ion 
constituent. 

Formal device denoting the molal 
activity coefficient of an ion 
constituent; eq (19). 

Thermodynamic composition co- 
efficient, eq (18). 

Stoichiometric coefficient of a 
neutral component in an elec- 
trode reaction, eq (12). 

Formal device denoting the 
chemical potential of an ion 
constituent. 

Difference in the chemical 
potential of an ion constituent 
between a and /5, eq (13). 

The chemical potential of a com- 
ponent represented in {m*}. 

The stoichiometry matrix, eq (4). 

The integral stoichiometric trans- 
ference number of a constituent 
in a cell fixed frame of reference ; 
eq (69). 

(Undefined) the absolute electro- 
static potential. 
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Equation (13). 
Conventionally 
section 9. 
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International practical temperature scale of 1948, text revision 
of 1960, H. F. Stimson, NBS Mono. 37 {I960) 10 cents. 
The International Practical Temperature Scale of 1948 is a 
text revision of the International Temperature Scale of 1948, 
the numerical values of temperatures remaining the same. 
The adjective ''Practical" was added to the name by the 
International Committee on Weights and Measures. The 
scale continues to be based upon six fixed and reproducible 
equilibrium temperatures to which values have been assigned, 
and upon the same interpolation formulas relating tempera- 
tures to the indications of specified measuring instruments. 
Some changes have been made in the text to make the scale 
more reproducible than its predecessor. The triple point of 
water, with the value 0.01 °(.' replaces the former ice point as 
a defining fixed point of the scale. It is also recommended 
that the zinc point, with the value 419.505 °C, be used in- 
stead of the sulfur j^oint. The recommendations include new 
information that ha^ become available since 1948. 

An evaluation of Kacser's second order Born approximation 
to the bremsstrahlung differential cross section, G. 8. OfeU, 
NHS Tech. Note 81 {PB161582) {1061) 75 cents. 
The second order term, as derived by C. Kacser [Proc. Roy. 
Soc. A253 (1959)1, of the Born approximation series for the 
bremsstrahlung cross section differential with respect to 
photon energy, photon direction, and final electron direction 
averaged over initial and summed over final polarization 
states has been put in a form that admits to numerical evalua- 
tion for the coplanar case. The results are valid for relati- 
vistic as well as non-relativistic incident electrons. 
Completely general (non-coplanar) expressions are included 
for the first order (Bethe-Heitler) and second order terms. 
However, the latter lias not been given in a form admissible 
to munerical evaluation. 

Tabulated values are given for the case of the incident elec- 
tron kinc^tic energy = 500 kev, plioton energy = 450 kev, and 
photon chrection 20° from the incident electron direction. 
For this case the second order term increases the differential 
cross section by about 50 percent for Al and about 200 per- 
cent for Au for the various angles of electron emission. A 
rough approximation to the integrated (over final electron 
direction) cross section is included and a comparison is made 
with experimental values. The inclusion of the second order 
term brings the theoretical cross section closer to the experi- 
mental values. 

A compilation of the physical equilibria and related properties 
of the hydrogen-carbon monoxide system, 1). 1]. Drayer and 
T. M. Flynn, NBS Tech. Note 108 {PB16160IJ) {1061) $2.25. 
Literature data have been used to calculate K factors (K = 
the ratio of vapor phase to condensed phase compositions) 
for the hvdrogen-carbon monoxide system over the range of 
68.2 to 122. 2°K and 10 to 225 atmospheres. K factors are 
presented graphically for eight isotherms over this range. 
PublislKKl data on the solid-vapor region are presented 
separately as composition versus pressure at constant temper- 
ature. 

A bibliography of approximately 450 references is also pre- 
sented on related properties for this system and for the pure 
components. 

A compilation of the physical equilibria and related properties 
of the hydrogen-helium system, 1). E. Drayer and T. M. 
Flynn, NBS Tech. Note 109 {PB161610) {1961) $1.25. 
Published data have been used to calculate K factors (K = 
the ratio of vapor phase to condensed phase compositions) 
for the helium-hydrogen system over the range of 17.4° to 



21.8°K and 2 to 32 atmospheres pressure. K factors are 
presented graphically for three isotherms over this range. 
A bibliography of approximately 290 references is also pre- 
sented on related properties for this system and for the pure 
components. 

A compilation of the physical equilibria and related properties 
of the hydrogen-nitrogen system, 1). 1]. Draver and T. M. 
Flynn, NBS Tech. Note 110 {PB161611) {1961) $1.75. 
Published data have been used to calculate K factors (K= 
the ratio of vapor phase to condensed phase compositions) 
for the hvdrogen-nitrogen system over the licjuid-vapor range 
of ()8.2 to 122. 2°K, and 10 to 225 atmospheres. K factors 
are presented graphically for eleven isotherms within this 
range. 

Published data on the solid-gas and solid-licpiid regions are 
presented separately as composition versus pressure at 
constant temperature. 

A l)il)liograj)liy of 250 references pertaining to the hydrogen- 
nitrogen system is included. 

The structure of B-type vibrational-rotational bands of an 
asymmetric rotor, H. C. Allen, Jr., Phil. Trans. Hoy. Soc. of 
London, Series A. Mathematical and Physical Sciences, No. 
1030, 253, 335 {1961). 

The structure of B-type vibrational-rotational ))ands of an 
asymmetric rotor is related to the band types of the two 
limiting symmetric rotors. The effect of the degree of 
asymmetry on the band structure is shown. The type of 
information which can be obtained from the analysis of a 
B-type band is pointed out. 

Complex conductivity of some plasmas and semiconductors, 

P. H. Fang, ./. Apvl- Pesearch Soc. B {The Hague) 9,51(1960). 
The complex conductivities of plasmas and semiconductors 
have been calculated for several cases where collision fre- 
quency can be expressed as a power function of the energy. 
From the result, some characteristic parameters of th(^ 
plasma originally investigated by Spitzer are estimated. 
The problem of determining the relaxation time from a non- 
symmetrical dispersion is discussed. 

Electric field distribution in a dense plasma, J. L. Jackson, 
Phys. Fluids 3, 927-93 1 {I960). 

A calculation is presented of the probability distribution 
function of the electric field at the center of an ion or atom 
in a plasma in the high ion density limit. In this limit, it 
is possible to take into account rigorously the effect of the 
Coulomb interactions on the distribution function. The 
distribution function in the high density limit is Gaussian. 
The Coulomb interaction decreases the mean square electric 
field by the multiplicative factor {l-\-kao-{- {kao)%)~\ 
where k is the reciprocal Debye length and ao the radius of 
the ion or atom at whose center the field is evaluated. 

Absorption spectrum and magnetic properties of osmium 
hexafluoride, J. C. Eisenstein, ./. Cheni. Phys. 34, 310-318 

{1961). 

This paper is devoted to a discussion of the optical absorp- 
tion spectrum and the susceptibility of CsFe. It is assumed 
that two 5d electrons, which move in a ligand field of octa- 
hedral symmetry, are responsible for the magnetic and 
optical properties of the molecule. The calculated positions 
of the energy levels depend on the strength of the ligand 
field, the Coulomb integrals and the spin-orbit coupling 
constant. These quantities are treated as parameters. 
If they are given appropriate values excellent agreement 
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with observed energy level positions is obtained. The 
calculated susceptibility is also in good agreement with 
experimental values. In order to fit the magnetic data it is 
necessary to assume that the orbital reduction factor is 
approximately 0.7. 

Effect of departures from local thermodynamic equilibrium 
on inferences on stellar atmospheric temperatures, R. N. 
Thomas, Book, Optical Spectrometric Measurements of High 
Temperatures, P. J. Dickerman, ed.y 14-26 {University of 
Chicago Press, Chicago, III,, 1961). 

A summary, representing the astrophysical approach to the 
spectrometric diagnosis of a gaseous atmosphere, for the 
Chicago conference on Optical Spectrometric Diagnostics of a 
High-temperature Gas. Problems of departure from Local 
Thermodynamic Equilibrium, in their bearing on the assign- 
ment of a local value of temperature in a gaseous atmosphere , 
are discussed. 

Thermal degradation of organic polymers, S. L. Madorsky, 
Polymer Set. 17, No. 7, (July 1961). 

When pyrolyzed in a vacuum or in a neutral atmosphere at 
about 200° to 1200° C, simple vinyl polymer chains break into 
molecular fragments of various sizes. The mechanism of the 
chemical reactions involved and the nature and relative 
amounts of the resulting products depend on the molecular 
structure of the polymer undergoing pyrolysis. Some 
polymers decompose into a molecular spectrum of fragments 
containing from 1 to 50 or more carbon chain atoms; some 
yield monomers exclusively; some yield a mixture of mono- 
mers and various other chain fragments; yet another group 
of polymers yield fragments which are unrelated to the 
structural unit of the chain. Generally, the higher the 
temperature and pressure of pyrolysis the greater the frag- 
mentation of the decomposition products. 

Infrared spectra of solid hydrocarbons at very low tempera- 
tures, J. J. Comeford and J. H. Gould, /. Mole. Spectroscopy 
5, 474 {I960). 

Mixtures of various hydrocarbons with argon were condensed 
on a Csl window cooled with liquid helium or hydrogen and 
the absorption spectra recorded in the region 700-3500 cm~i. 

Microwave Zeeman effect of free hydroxyl radicals: ^\l\i2 
levels, H. E. Radford, Phys. Rev. n%, 114 {1961). 

Two new paramagnetic resonance spectra of the free O^^H 
radical have been observed at 3 cm wavelength in the products 
of an electric discharge in water vapor. The spectra arise from 
x\-type doubling transitions in the ^iji^ J = 3/2 and '^llh, 
J =5/2 levels. Molecular magnetic moments and hyperfine 
structure constants measured from the spectra lead to refine- 
ments in the theory developed earlier to account for similar 
measurements on '^11^/2 levels. A consistent analysis of the 
combined ^H^ and 2II3/2 data is then possible, and yields a 
much better understanding of the molecular Zeeman effect 
and hyperfine structure. 

Kinetic isotope effects in the reaction of methyl radicals 
with ethane, -do, and ethane-1, 1, \-dz, J. R. McNesbv, 
J. Phys. Chem. 64, No. 11, 1671 {Nov. 1960). 

Photolyses of acetone-Je in the presence of a mixture of C2H6 
and C2D6 and in the presence of CH3CD3 were carried out to 
measure kinetic isotope effects for the reactions 

(1) CD3+C2H6 -^CDsH + CsHs -> 

(2) CDg+CHsCDs-^CDsH + CHsCDs-^ 

(3) CD3 + CH3CD3->CD4 + CH3CD2 

(4) CD3+C2D6 -^CD4+C2D5 

The conclusion drawn from the results is that the relative rate 
constants are within experimental error: 

ki/h = exp{ldOOlRT) 

k2/h = exp{ldOO/RT) 

ki/k2 = 2.0 

Theory of vibrational relaxation in liquids, R. Zwanzig, 
J. Chem. Phys. 34, 1931-1935 {1961). 

A new formulation of the theory of vibrational relaxation, 
based on Zener's semi-classical approximation, is presented 
here. The relaxation rate is shown to be proportional to the 



spectral density of the force exerted on the oscillator by its envi- 
ronment. The isolated binary collision theory is derived, but 
only with the condition that the collision frequency is much 
smaller than the oscillator frequency. This requirement is 
not satisfied in a liquid: we conclude that Litovitz's appli- 
cation of the isolated binary collision theory to liquid CS2 is 
not justified. A possible relation between vibrational relax- 
ation and the self-diffusion coefficient in a liquid is discussed. 

Nomenclature of carbohydrates, R. S. Tipson, J. Chem. 
Documentation 3, No. 8, 3-7 (1961). 

In naming a carbohydrate, the name chosen should not 
conflict with the general principles of established organic 
nomenclature; it should concisely show the number of carbon 
atoms, the number and configuration of all asymmetric 
carbon atoms, the functional groups, the size of any sugar- 
ring, and the anomeric form_, if any. 

The present paper describes methods which have evolved for 
naming sugars having either one or two aldehydic or ketonic 
functions, or an aldehydic plus a ketonic function. The 
system is then extended to derivatives having either one or 
two carboxyl groups, an aldehydic plus a carboxyl function, 
or a ketonic plus a carboxyl function. Naming of derivatives 
of such compounds are also discussed. 

The oxide films formed on copper single crystal surfaces in 
water 11. Rate of growth at room temperature, J. Kruger, 
/. Electrochem. Soc. 108, No. 6, 503 {June 1961). 
Measurements utilizing polarized light were made of the in- 
crease in film thickness with time on {100}, jllO), {111}, 
and {311} surfaces of a copper single crystal immersed in 
water in equilibrium with atmospheres containing either pure 
oxygen or oxygen-helium mixtures containing 1%, 10%, or 
20% oxygen. With the pure oxygen atmosphere a CU2O 
film rapidly reached a limiting thickness that remained 
constant for 1.5 to 2 hours after which time the formation of 
CuO was observed. Stirring increased the time required to 
observe the formation of CuO. Similar behavior was ob- 
served for 10% and 20% oxygen atmospheres, the limiting 
film thickness of CU2O and the time required to observe the 
formation of CuO increased with decreasing oxygen content 
in the atmosphere. When a 1% oxygen atmosphere was 
used, parabolic film growth was observed. In all cases the 
{111} plane exhibited the lowest growth rate. 

Line shape of ultraviolet absorption in solid noble gases, 

P. H. E. Meijer, J. Chem. Phys. 34, 2078-2082 {1961). 
The absorption of ultraviolet in solid noble gases shows lines 
that are a few inverse centimeters broad. An attempt is 
made to explain and calculate the width and shape of these 
lines on the basis of a tentative energy versus configuration, 
curve. It is shown that the width has the proper order of 
magnitude. The line shape can be easily calculated at the 
long wave length end. It should show an exponential tail. 
If we inspect the expression for small values of the energy, 
we find a rapidly rising exponential, which would give rise to 
a smeared out edge appearance, provided the normalization 
constant (which cannot be determined in the approximation 
used) is relatively insensitive to variations in the energy. 
The line will have a nonsymmetrical shape. 

Exponential temperature dependence of Young's modulus 
for several oxides, J. B. Watchman, Jr., W. E. Tefft, D. G. 
Lam, Jr., and C. S. Apstein, Phys. Rev. 123, 1754 {June 1961). 
Young's modulus was measured over the temperature range 
77 °K to 850 °K by an accurate resonance technique. Data 
are presented for single crystals of aluminum oxide with 
various orientations of the crystallographic axes and for 
polycrystalline aluminum oxide, thorium oxide, and mag- 
nesium oxide. The results show that the range of validity 
of a T^ temperature dependence predicted by theory must be 
quite small. The temperature dependence is very well de- 
scribed over the whole temperature range by Texp(— To/^) 
where To is an empirical parameter. 

The X- or gamma-ray energy absorption of transfer co- 
efficient: Tabulations and discussion, R. T. Berger, Radiation 
Research 15, No. 1, 1-29 {July 1961). 

The energy absorption coeflficient is defined and its signifi- 
cance for the calculation of absorbed dose is discussed. 



98 



Values of the energy absorption coefficient are tabulated for 
X-ray photon energies between 3 kev and 10 Mev, for various 
materials. This tabulation is an extension of previous work 
of CJ. R. White Grodstein and takes into account the latest 
values of the relevant cross sections. Auxiliary data are 
given to facilitate adaptation of the energy absorption 
coefficient to practical applications. The uncertainties of the 
ta})i dated values are discussed. 

Adsorption of nitrogen and argon on mineralogical graphite 
and diamond at 77 and 90° K, J. de Dios Lopez-Gonzalez, 
F. (J. Garpenter and V. R. Deitz, J. Phys. Chem. 65, No. 7, 
1112-19 (July 1961). 

Adsorption isotherms of argon and nitrogen on mineralogical 
graphite and diamond powder were determined in the relative 
pressure range from 10"^ to about 1. The sample was out- 
gassed at 450° C for 18 hours before each of a series of measure- 
ments. Upon the first exposure to the adsorbate, there was in 
many cases a very slow adsorption. This was followed by 
a much higher rate of adsorption and finally by attainment of 
a steady state. An induction period in physical adsorption 
is postulated to explain this phenomenon. The isotherms on 
diamond followed the Freundlich equation over a very wide 
I^ressure range and the B.E.T. relationship from 0.01 to 0.4 
relative pressure. The isotherms on graphite showed a 
definite step at a relative pressure range betw^een 10~^ and 
lO""*, and indicated another step at about 10"^ The iso- 
therms on graphite followed the B.E.T. relationship in the 
relative pressure range 0.001 to 0.1. 

Absorption of carbon dioxide by solutions of 2-amino-2- 
(hydroxymethyl)-l, 3-propanediol, R. G. Bates and 11. B. 
Hetz(;r, Anal. Chem. 33, No. U, 1285 {Aug. 1961). 
Tr?'s(hydroxy methyl) aminomethane, or 2-amino-2-l, 3-pro- 
panediol, is finding wide use in acidimetry and in pH control. 
This useful base and its aqueous solutions have been reported 
to be stable in air. Nevertheless, an appreciable change of 
pH can occur if these solutions are exposed to atmospheres 
containing normal amounts of carbon dioxide. 

Quantum theory of interference effects in the mixing of light 
from phase-independent sources, U. Fano, Am. J. Phys. 29, 

No. 8, 539-545 {August 1961). 

Correlations in space and time of the light intensity from 
extended sources, observed by Brown and Twiss and by 
Forrester, were analyzed by these authors primarily in terms 
of field oscillations. This paper discusses an atomic process 
which brings out the intensity correlations, namely, the 
photoionization of a pair of atoms follow^ing photon emission 
by another pair of independently excited atoms. The cal- 
culated probability of this process depends sinusoidally on 
the relative positions of the four atoms and also, when the 
source atoms emit different frequencies, on the time interval 
between the photoionizations, in agreement with the macro- 
scopic treatment. The oscillations arise from an interference 
of probability amplitudes which is not affected by the random 
phases of the source atoms. The calculation follows a stand- 
ard approach but involves some novel detail. 

Heats and volumes of mixing in several C12 hydrocarbon 
systems, R. S. Jessup and C. L. Stanlev, /. Chem. Eng. Data 
6, No. S, 368-371 {July 1961). 

Measurements have been made of the heats and volumes of 
mixing of two binary systems, one of which was composed of 
two kerosines, and the other of the C12 hydrocarbons bicyclo- 
hexyl and n-dodecane. In addition, measurements of volume 
of mixing were made on four other binary systems, in each of 
which one component was n-dodecane while the other was, 
in turn, n-heptylcyclopentane, n-hexylcyclohexane, n-hexyl- 
benzeno, and phenylcyclohexane. 

Low-angle X-ray diffraction of crystalline nonoriented poly- 
ethylene and its relation to crystallization mechanisms, 

L. Mandelkern, A. S. Posner, A. F. Diorio, and D. E. Roberts, 
/. Appl. Phys. 33, No. 8, 1509-1517 (Aug. 1961), 

X-ray diffraction maxima at low angles have been observed in 
crystalline but nonoriented linear polyethylene, the crystalli- 
zation process being conducted from the melt of the pure 
undiluted polymer. Several orders of diffraction are ob- 
served in favorable cases, and the spacings corresponding to 



the first-order reflections range from 150-850 A. The values 
of the maxima depend on the mode of crystallization. Previ- 
ous assertions that in such systems the maxima are limited to 
100-200 A are shown to be very restrictive and typical only 
of crystallization processes conducted at very large under- 
cooling. Major attention is focused on the properties of 
specimens crystallized isothermally at relatively low values of 
the undercooling. The spacings are very sensitive to the 
crystallization temperature in this range. The highest 
values are observed at low undercooling, and substantial 
decreases occur as the temperature^ is lowered. Goncomi- 
tantly, the density observed after isothermal crystallization 
significantly decreases with a lowering of the crystallization 
temperature. 

The fact that a periodicity can be developed in such systems, 
the magnitude of the maxima, and their dependence on the 
crystallization temperature, is explicable by the application of 
nucleation theory. It is assumed that subsequent to the 
formation of critical-size nuclei from a bundle of polymer 
chains, crystal growth along the chain direction is severely 
retarded, while in the transverse direction essentially unim- 
peded crystallization occurs. From the observed tempera- 
ture coefficient of the low-angle spacings, the ratio of the 
excess free energy (resulting from the junction of crystalline 
and amorphous regions at the crystallite ends) to the bulk 
enthalpy of fusion is found to be 2.6. The magnitude of this 
ratio receives confirmation from another type of experiment. 

Density fluctuations and heat conduction in a pure liquid, 

R. E. Nettleton, Phijs. Fluids 4, No. 1, 74 {Jan. 1961). 
Two components of the heat flux in a pure liquid are identi- 
fied — one carried by sound waves and the other by diffusing 
molecules. Coupled relaxation equations are obtained for 
these components, both depending on the density gradient 
Vp, and this dependence is interpreted thermodynamically 
as a further coupling of these equations with the equation 
determining {d/bt) (Vp). The latter can be immediately 
written down with the aid of Onsager's theorem and an 
inertial principle developed in a previous paper. The thermo- 
dynamic interpretation of these rate equations also leads to 
an explicit expression for the coefficient of (Vp)^ in a Taylor 
expansion, about equilibrium, of the Helmholtz free energy. 
This result is compared with similar free energy terms assumed 
to exist by other authors. 

Low temperature properties of cerous magnesium nitrate, 

R. \\ Hudson, R. S. Kaeser, and H. E. Radford, Proc. Vllth 
Intern. Conf. Low Temperature Physics, University of Toronto, 
Canada, p. 45-48 {1950). 

A variety of measurements have been made on cerous magne- 
sium nitrate between room temperature and the region of 
millidegrees Kelvin. In the range 10"^ to 1° K the entropy 
temperature relation varies from specimen to specimen and 
at the high temperature end the curve for any one sample 
does not entropolate to coincidence with that derived from 
spin-lattice relaxation measurements in the l°-4° K region. 
The spectroscopic splitting factor is probably exactly zero 
along the c-axis. Evidence for the presence of low-lying 
doublets is supplied by departures of the X-T relation from 
the Curie Law in the 4°-300° K range. 

Other NBS Publications 

Journal of Research 65B (Math, and Math. Phys.) No. 4 
(Oct.-Dec. 1961) 75 cents. 

Physical entities and mathematical representation. G. H. 

Page. 
On the range of a fleet of aircraft. A. J. Goldman. 
Measurement of wave fronts without a reference standard: 

Part I. The wave-front-shearing interferometer. J. B 

Saunders. 

On the evaluation of the function 'i>(X) =77—. I e"'°"+^"d2^ 

ZtTZ Jff-ico 

for real values of X. W. Borsch-Supan. 
Analyticity and probability properties of one-dimensional 
Brownian motion. A. Ghaffari. 
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Some higher order integral identities with appUcation to 

bounding techniques. J. H. Bramble and B. E. Hubbard. 
A priori bounds in the fir«t boundary vahie problem in 

elasticity. J. H. Brambk\ and L. E. Payne. 
Journal of Research 66D (Radio Prop.) No. 1 (Jan.-Feb. 
1962) 70 cents. 
A survey of very wide band antennas — 1945 to the present 

J. D. Dyson. 
Numerical investigation of the equivalent impedance of a 

wire grid parallel to the interface between two media. 

T. Larsen. 
Current on and input impedance of a cylindrical antenna. 

Y. M. Chen and J. B. Keller. 
Kadar corner reflectors for linear or circular polarization. G. 

Latmiral and A. Sposito. 
On the theory of wave propagation through a concentrically 

stratified troposphere with a smooth profile. H. Bremmer. 
On the propagation of VLF and ELF radio waves when the 

ionosphere is not sharply bounded. J. R. Wait. 
Fields of electric dipoles in sea water — the earth-atmosphere- 
ionosphere problem. W. L. Anderson. 
Reflection of electromagnetic waves from thin ionized gaseous 

layers. F. H. Northover. 
Reflection and transmission of radio waves at a continuously 

stratified plasma with arbitrary magnetic induction. J. R. 

Johler and J. D. Harper, Jr. 
On the diffraction of spherical radio waves by a finitely 

conducting spherical earth. L. C. Walters and J. R. 

Johler. 
An approximate full wave solution for low frequency electro- 
magnetic waves in an unbounded magneto-ionic medium. 

W. C. Hoffman. 
VHF radio propagation data for the Cedar Rapids-Sterling, 

Anchorage-Barrow, and Fargo-Churchill test paths, April 

1951 through June 1958. G. R. Sugar and K. W. Sullivan. 



Stopping powers for use with cavitv chambers, NBS Handb. 

79 (1961), 35 cents. 
Hydraulic research in the United States, H. K. Middleton, 

NBS Misc. Publ. 238 (1961), $1.25. 
Fractional factorial designs for experiments with factors at 

two and three levels, W. S. Connor and S. Young, NBS 

Applied Math. Series 58 (1961), 40 cents. 
Quarterlv radio noise data, March, April, Mav 1961, W. Q. 

Crichlbw, R. T. Disney, and M. A. Jenkins, NBS Tech. 

Note 18-10 (PB151377~10) (1961) $1.50. 
Mean electron densitv variations of the quiet ionosphere. 

No. 5— Julv 1959, ^ J. W. Wright, L. R. Wescott, and 

D. J. Brown, NBS Tech. Note 40-5 (PB151399-5) (1961) 

$1.50. 
Flux switching mechanisms in ferrite cores and their depend- 
ence on core geometrv, G. W. Reimherr, NBS Tech. 

Note 90 (PB161591) (1961) $1.25. 
An experimental studv of beta decav using the radiations from 

oriented nuclei, D. D. Hoppes, NBS Tech. Note 93 

(PB161594) (1961) $1.50. 
Bibliography on meteoric radio wave propagation, W. Nupen, 

NBS Tech. Note 94 (PB161595) (1961) $2.75. 
The integrated starlight over the skv, L. R. Megill and F. E. 

Roach, NBS Tech. Note 106 (PB161607) (1961) $2.00. 
Load carrying capacity of gas-lubricated bearings with 

inherent orifice compensation using nitrogen and helium 

gas, H. Sixsmith, W. A. Wilson, and B. W. Birmingham, 

NBS Tech. Note 115 (PB161616) (1961) $1.00. 
Matrix elements in the forbidden beta decay of Ce^^\ D. D. 

Hoppes, E. Ambler, R. W. Hay ward, and R. S. Kaeser, 
^ Phys. Rev. Letters 6, 115-118 (1961). 
The electromagnetic fields of horizontal dipole in the presence 

of a conducting half space, J. R. Wait, Can. J. Phys. 39, 

No. 7, 1017-1028 (July 1961). 
Simple adiabatic demagnetization apparatus, V. D. Arp and 

R. H. Kropschot, Rev. Sci. Instr., Note 33, 217-218 (Feb. 

1961). 
4-Nitrophenyl esters, carbonates, and bi-carbonates as leather 

fungicides, S. Dahl and A. M. Kaplan, J. Am. Leather 

Chemists Assoc. 55, 480 (Sept. 1960). 
The organization of international intercomparisons of radio- 
activity standards, with special reference to such measure- 



ments of NBS standards, L. M. Cavallo and W. B. Mann, 
Proc. 1959 Symp. Metrology Radioactive, IAEA, Vienna, 
^ 117 (1960). 

The physical environment as affected by radiation, D. M. 
Gates, Proceedings of the 50th Anniversary Celebration of 
the University of Michigan Biological Station (1960). 

Second revision American Dental Association specifications No. 
4 for dental inlay casting wax, J. W. Stanford, K. V. Weigel, 
and J. W. Paffenbarger, J. Am. Dental Assoc. 63, 45 (Jan. 
1961). 

Measurement of contrast in the aerial image, F. W. Rosberrv, 
Photogrammetric Eng. 155-159 (March 1961). 

An interferometric instrument for the rapid measurement of 
small diameters, D. H. Blackburn, Rev. Sci. Instr. 33, No. 2, 
137 (Feb. 1961). 

Vapor-phase photolvsis of formic acid, R. Gordon, Jr., and 
P. Ausloos, J. Chem. Phys. 65, 1033 (1961). 

Binary silicate glasses in the study of alkali-aggregate reaction, 
R. G. Pike, D. Hubbard, and E. S, Newman, Highway Re- 
search Board BuU. No. 275, 39-44 (1960). 

Phvsical measurement — challenge to science and engineering, 
A. V. Astin, SPE J. 17, No. 5. 455-458 (May 1961). 

Low-temperature thermocouples I. Gold-cobalt or constantan 
versus copper or '^normal" silver, R. L. Powell, M. D. 
Bunch, and R. J. Corruccini, Cryogenics 1, 139-150 (Mar. 
1961). 

Dielectric properties of polvamides, A. J. Curtis, J. Chem. 
Phys. 34, No. 5, 1849-1850 (May 1961). 

Natural and svnthetic rubbers: Review of test methods, F. J. 
Linnig, M. Trvon, and E. J. Parks, Anal. Chem. 33, 1272 
(Apr. 1961). 

Simple calculus for all-dielectric interference filters, K. D. 
Mielenz, J. Opt. Soc. Am. 50, No. 10, 1014 (Oct. 1960). 

Magnetic susceptibility of tetratonal titanium dioxide, F. E. 
Senftle, T. Pankey, and F. A. Grant, Phys. Rev. 130, 820 
^ (Nov. 1960). 

Equipment and procedures for the evaluation of total hemi- 
spherical emittance, J. C. Richmond and W. N. Harrison, 
Bull. Am. Coram. Soc. 39, No. 11, 668 (Nov. 1960). 

Limitations of radiosonde punch-card records for radio-mete- 
orological studies, B. R. Bean and B. A. Cahoon, J. 
Geophys. Research 66, No. 1, 328-331 (Jan. 1961). 

Evaluations of some oxide glasses for use as infrared materials, 
G. W. Cleek and T. G. Scuderi, Proc. Infrared Information 
Symp. 5, No. 4 (Aug. 1960). 

Improvement of some of the properties of sole leathers by im- 
pregnation with polvmers, J. R. Kanagv, J. Am. Leather 
Chemists Assoc. LVI, No. 7, 322-342 (July 1961). 

A transistor frequencv meter, F. R. Bretemps and S. Saito, 
Electronic Inds. 19,^ No. 10, 196-198 (Oct. 1960). 

Studies on the fiame-spraving of aluminum oxide, D. G. 
Moore, Metal Finishing J. 6, No. 70, 397 (Oct. 1960). 

Self-qualification of laboratories, A. T. McPherson, Instr. and 
Control Systems 34, No. 7. 1265 (July 1961). 

Four-color achromats and superchromats, R. E. Stephens, J. 
Opt. Soc. Am. 50, No. 10, 1016-1019 (Oct. 1960). 

Bounds for determinants with positive diagonals, E. V. 
Haynsworth, Trans. Am. Math. Soc. 96, 395-399 (1960). 

Refining measurements by capacitance techniques, F. K. 
Harris and R. D. Cutkosky, ISA — J. Instr. Soc. Am. 8, 
No. 2, 63-66 (Feb. 1961). 

Standards for plastic pipe fittings, F. W. Reinhart, SPE J. 17, 
160 (Feb. 1961). 

Post office mechanization, B. M. Levin, M. C. Stark, and P. C. 
Tosini, Elec. Engr. 80, 105-110 (Feb. 1961). 

Comparison of observed tropospheric refraction with values 
computed from the surface refractivity, B. R. Bean, IRE 
Trans, Ant. Prop. AP-9, 415-416 (July 1961). 

The precision measurement of transformer ratios, R. D. 
Cutkosky and J. Q. Shields, IRE Trans, Instrumentation 
1-9, No.^2, 243-250 (Sept. 1960). 

Studies of environmental factors in a family-size underground 
shelter, P. R. Achenbach, F. J. J. Drapeau, and C. W. 
Phillips, Report OCDM-NBS-60-1 issued Office of Civil 
Defense MobiUzation (Mar. 1961). 

Electrodes for pH measurements, R. G. Bates, J. Electro- 
analytical Chem. 3, 93 (1961). 

Phase-shift method for one-dimensional scattering, A. H. 
Kahn, Am. J. Phys. 39, 77 (1961). 
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The classical field theories, C. Truesdell and R. Toupin, 
Handb. der Physik, Springer, Berlin, III, No. 1, 226-793 
(1960). 

A precision RF attenuation calibration system, C. M. Allred 
and C. C. Cook, IRE Trans. Instrumentation 1-9, No. 2, 
268-274 (Sept. 1960). 

A vibniting sample magnetometer, N. V. Frederick, IRE 
Trans. Instrumentation 1-9, No. 2, 194-196 (Sept. 1960). 

Standards development by ASTM committee D-20 on 
plastics, F. Reinhart, SPE J. 16, 800 (Aug. 1960). 

A nuclear magnetic resonance study of sny-onti isomerism in 
kctoxinies, E. Lustig, J. Phys. Chem. 65, 491 (1961). 

A standard current transformer and comparison method — 
A basis for establishing ratios of currents at audio frequen- 
cies, H. L. Dunfee, IRE Trans. Instrumentation 1-9, No. 2, 
231-236 (Sept. 1960). 

Evaluating the freezing-and-thawing durability of concrete 
by laboratory tests in the U.S.A., H. T. Arni and R. L. 
Blain(^, Preliminary Report on Symposium on Durability 
of Concrete of RILEM, Prague, 'Czechoslovakia, 218-231 
(July 1961). 

Spectroscopic standard samples of titanium and high- 
temperature alloys, R. E. Michaelis, J. Appl. Spectroscopy 
15, 7 (1961). 

Analytical laboratory problems, J. Mandel, Trans. Metro- 
politan Conf. Am. Soc. Quality Control, 99-107 (Feb. 
1960). 

Some exporinKmtal aspects of nuclear orientation, E. Ambler, 
Proc. lOth Inter. Congress of Refrigeration, Copenhagen, 
Denmark, 1959, Progress in Refregeration in Science and 
Technology I, 195-198 (Pergamon Press, London, England, 
1960). 

Rapid method for interpolating refractive index measure- 
ments, O. N. Stavroudis and L. E. Sutton, J. Opt. Soc. 
Am. 51, No. 3, 368-370 (Mar. 1961). 

Comparison of four different methods of determining drying 
shrinkage of concrete masonry units, J. O. Bryson and 
D. Watstein, J. Am. Concrete Inst. 58, No. 2, 163-184 
(Aug. 1961). 

The retannage of chrome-tanned leather with vegetable 
tannins, J. R. Kanagy, J. Am. Leather Chemists Assoc. 
LVI, No. 6, 252-272 (June 1961). 

Improved electron filter lens, J. A. Simpson and L. Marton, 
Rev. Sci. Instr. 3^, No. 7, 802-803 (July 1961). 

Static and dynamic calibrations of pressure measuring in- 
struments at the National Bureau of Standards, E. C. 
Llovd and D. P. Johnson, Automatic and Remote Control 
Proc. 1st Intern. Congress of IF AC, 274-280 (1960). 

X-ray diffraction study of acellular teleost bone, M. L. Moss 
and A. S. Posner, Nature 188, 1037 (Dec. 1960). 

The NBS absolute gravity experiment, D. R. Tate, Proc. 2d 
Air Force Cambridge Research Center Military Geodesy 
Seininar, GRD Research Notes, 40 (Jan. 12-14, 'i960). 

Adsorption of polyesters and other polymers on glass and 
other substrates, R. R. Stromberg and G. M. Kline, 
Modern Plastics 38, 123 (Apr. 1961); 38, 241 (May 1961); 
Poliplastic 9, 15 (Jan.-Feb. 1961). 

Tunnel diode large-signal simulation study, S. B. Geller and 
P. A. Mantek, Proc. IRE 49, No. 4, 803 (Apr. 1961). 

Determination of carboxyl in cellulose — comparison of various 
methods, W. K. Wilson and J. Mandel, Tappi 44, No. 2, 
131 (1961). 

A recalibration of NBS carbon-14 standards by Geiger-Miiller 
and proportional gas counting, W. B. Mann, H. H. Seliger, 
W. F. Marlow, and R. W. Medlock, Rev. Sci. Instr. 31, 
No. 7, 690 (1960). 

Computer for weather data acquisition, P. Meissner, J. 
Cunningham, and C. Kettering, Proc. Eastern Joint Com- 
puter Conf. 18, No. 2.1, 57-66 (Dec. 1960). 

The solar spectrum from 2635 to 2085 A, H. H. Malitson, 
J. D. Purcell, R. Tousev, and C. E. Moore, Astrophvs. J. 
133, 746 (1960). 

Winter thermal radiation studies in Yellowstone Park, D. M. 
Gates, Science 134, 32-35 (July 7, 1961). 

Vapor phase 7-radiolysis of acetone, L. J. Stief and P. Ausloos, 
J. Phys. Chem. 65, 877 (1961). 

Amine buffers for pK control, R. G. Bates, Ann. N. Y. 
Acad. Sci. 92, 341 (1961). 

Cryogenic adhesive properties of bisphenol-A epoxy resins, 



M. J. Hiza and P. L. Barrick, SPE Trans. 1, 73-79 (A])r. 

1961). 
Resonance characteristics of a corrugated cylinder excited 

by a magnetic dipole, J. R. Wait and A. M. Conda, IRE 

Trans. Ant. Prop. AP-9, 330-333 (July 1961). 
Recent work with compensated internal gas counters for the 

standardization of gaseous radionuclides, W. B. Mann, 

Proc. Svmp. Metrology of Radioactive IAEA, Vienna, 

307 (1960). 
Vanadium oxytrichloride [VanadyUV) chloride], R. B. 

Johannesen, Book, Inorganic Syntheses, E. G. Rochow, 

Editor-in-Chief (McGraw-Hill Book Co., New York, N.Y.), 

IV. Chapt. 5B, No. 38, 119-120 (1960). 
Concerning the potential refractive index and the molecular 

refractivity, B. R. Bean and J. D. Horn, J. Meteorol. 18, 

No. 3, 427-428 (June 1961). 
Standard samples and related materials for spectrochemical 

analysis, R. E. Michaelis, Am. Soc. Testing Materials 

Spec. Tech. Publ. No. 58-D, 116 (1960). 
Aromatic fluorocarbons: Thermal stability and synthesis, 

L. A. Wall, Proc. 6th Joint Army-Navy-xAir Force Elastomer 

Conf., Oct. 18-20, 1960, I. 150 (U.S. Army Quartermaster 

Research & Engineering Command, Natick, Mass.). 
The ionization constant of hvdroxvlamine, R. A. Robinson 

and V. E. Bower, J. Phys. Chem'. 65, 1279-1280 (1961). 
The calculation of the energy levels of nearly symmetric 

rotors, S. C. Wait, Jr., J. Mol. Spectroscopy 6, 276 (1960). 
Multiple layer insulation for cryogenic applications, R. H. 

Kropschot, Cryogenics 1, 171-177 (Mar. 1961). 
Extinguishment effectiveness of some powdered materials on 

hydrocarbon fires, T. G. Lee and A. F. Robertson, Fire 

Research Al)stracts and Reviews, Natl. Acad. Sci.-Natl. 

Research Council Publ. 786 (June 1961). 
Forty years of dental research at the National Bureau of 

St'andards, W. T. Sweeney, Capital Chemist 10, 228 (Get. 

1960). 
A stainless steel for standard weights, S. J. Rosenberg and 

T. P. Royston, Materials Research & Standards (ASTM 

Bull.) 1, No. 1, 21 (Jan. 1961). 
Long-term working stress of thermoplastic pipe, F. W. 

Reinhart, SPE J., 751-754 (Aug. 1961). 
NBS atomic frequency standards, R. C. Mockler and R. E. 

Beehler, Proc. 14th Annual Symp. Frectuency Control, 

Atlantic City, N.J., May 31-June 3, 1960, p. 298-309 

(Signal Corps Eng. Laboratories, Fort Monmouth, N.J., 

1961). 
A precision RF power transfer standards, P. A. Hudson, 

IRE Trans. Instrumentation 1-9, 280-283 (Sept. 1960). 
The fluorescence and phosphorescence of trifluoroacetone 

vapor, P. Ausloos and E. Murad, J. Am. Chem. Soc. 83, 

1327 (1961). 
The morphology of sporadic E, E. K. Smith, Jr., Proc. 

Xlllth Gen. 'Assembly IJRSI (London, England, 1960), 

Natl. Acad. Sci.-Natl. Research Council Pul^l., Report No. 

880, 218 (1961). 
How to determine stack load limits, R. S. Wyly, Contractor, 

p. 19-20 (Aug. 15, 1961). 7-Irradiation of liquid and 

solid oxygen, D. W. Brown and L. A. Wall, Phys. Chem. 

65, 915 (1961). 
Testing of ball bearings with five difi'erent separator materials 

at 9200 RPM in hquid nitrogen, J. A. Brennan, W. A. 

Wilson, R. Radebaugh, and B. W. Birmingham, Am. Soc. 

Mech. Engrs., No. 61-LITBS-18, 1-8 (Apr. 1961). 
The CRPL electron density profile program: Some features 

and early results, J. W. Wright, Proc. URSI-AGI Symp. 

(Brussels, Belgium, 1959), Book, Some Ionospheric Results 

Obtained During IGY, Ed. W. J. G. Beyon (Elsevier 

Publ. Co., Princeton, N.J.), 215 (1960). 
Digest of a general description of D-C digital voltmeters, C. 

Stansbury, AIEE Trans. Paper 61-721 (Apr. 1961). 
Determination of wax in paper, W. L. Marrow, Tappi 44, 

No. 2, 120 (1961). 
Mortar for cavity walls, C. C. Fishburn, Building Research 

Institute, 1960 Spring Conf. Program on Insulated Cavity 

Walls, Natl. Acad. Sci.-Natl. Research Council Publ. 793, 

37 (1960). 
Variations of [01] 5577 A emission in the upper atmosphere, 

F. E. Roach, Ann. Geophys. 17, 172-180 (Apr.- June 1961). 
Continuously operating He^ refrigerator for producing 
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temperature down to }^° K, E. Ambler and R. B. Dove, 
Rev. Sci. Instr. 32, No. 6, 737-739 (June 1961). 

Location of the plane of best average definition for airplane 
camera lenses, F. E. Washer and W. P. Tayman, Photo- 
grammetric Eng. 26, No. 3, 475-488 (1960). 

A review of the hydrauhcs of circular sewers in accordance 
with the Manning formula, R. S. Wyly, Yearbook of Am. 
Soc. of Sanitary Engineering 38, 154 (1960). 

Non-reciprocity and time-reversal in microwave circuits, H. A. 
Fowler, Academie royale de Belgique (Classe Des Sciences) 
V-46, No. 11, 963-977 (Dec. 1960). 

A screening method for large information retrieval systems, 
R. T. Moore, Proc. Western Joint IRE-AIEE-ACM Com- 
puter Conf., Los Angeles, Calif. 6.3, 259-274 (1961). 

Humidity standards, A. Wexler, Tappi 44, No. 6, 180a (June 
1961). 

Tchebycheff approximations by ah''-\-c, J. R. Rice, J. Soc. Ind. 
AppL Math. 8, No. 4, 691-702 (1960). 

Ionospheric ''forward" scattering, D. K. Bailey, Proc. Xlllth 
Gen. Assembly of URSI (London, England, 1960), Natl. 
Acad. Sci.-Natl. Research Council Publ. Report No. 880, 
281 (1961). 

Apparatus for controlled slack quenching, N. L. Carwile, 
M. R. Meyerson, and S. J. Rosenberg, Materials Research 
and Standards (ASTM Bull.) 1, No. 7, 532-536 (July 1961). 

A high-resolution ammonia-maser-spectrum analyzer, J. A. 
Barnes and L. E. Heim, IRE Trans. Instrumentation I-IO, 
4-8 (June 1961). 

Electron spin resonance spectra of free-radical intermediates 
formed by reaction of polystyrene with atoms of hydrogen 
and deuterium, R. B. Ingalls and L. A. Wall, J. Chem. 
Phys. 35, No. 1, 370-371 (July 1961), 

Reinhart, F. M., Tw^enty-year atmospheric corrosion investi- 
gation of zinc-coated and uncoated wire and wire products, 
Am. Soc. Testing Materials Spec. Tech. Publ. 290, 1-141 
(June 1961). 

Time standards, A. G. McNish, Instr. and Control Systems 
33, 1340-1341 (Aug. 1960). 



Stress-strain relationships in yarns subjected to rapid impact 
loading; Part 7. Stress-strain curves and breaking-energy 
data for textile yarns, J. C. Smith, P. J. Shouse, J. M. 
Blandford, and K. M. Towne, Textile Research J. 31, No. 8, 
721-734 (Aug. 1961). 

V . Mutual interference between surfaces and satellite communi- 
cation systems, B. Hartman, M. T. Decker, eh. 5, p. 111- 
157, Book, Frequency Allocations for Space Communica- 
tions (Report of the Joint Tech. Adv. Comm. of IRE-EIA 
to the Federal Communications Commission, Mar. 1961). 

Some observations of copper deposits on single crystals of 
copper, I. Giron and F. Ogburn, J. Electrochem. Soc. 108, 
No. 9, 842-846 (Sept. 1961). 

Binary oxide systems, H. R. Hoekstra and H. S. Parker, ch. 6, 
sec. 6, p. 271-304, Book, UO2, Properties and Nuclear 
Applications (U. S. Superintendent of Documents, Govern- 
ment Printing Office, Washington 25, D. C, Aug. 1961). 

Wind resistance of asphalt shingle roofing, W. C. Cullen, 
Natl. Acad. Sci.-Natl. Research Council, 33-42 (1961). 

Pulse voltage comparator measures height of positive or nega- 
tive pulses, 0. B. Laug, Electronics 34, No. 36, 70-71 
(Sept. 1961). 

Spectroscopic observations of the stratosphere and mesosphere 
D. M. Gates, Instr. Soc. Am. J. Conf. Preprint No. 2-NY60 
(Sept. 1960). 

Supercoohng in the barnyard, C.^P. Saylor, The Capital 
Chemist 10, No. 9, 292-293 (Dec. 1960). 

References to contemporary papers on acoustics, W. Koidan, 
Acoust. Soc. Am. 33, No. 5, 705-716 (May 1961). 

'^Publications for which a price is indicated {except for Tech- 
nical Notes) are available only from the Superintendent of Docu- 
ments, U. S. Government Printing Office, Washington 25, D. C. 
(foreign postage, one-fourth additional). Technical Notes are 
available only from the Office of Technical Services, U. S. Depart- 
ment of Commerce, Washington 25, D. C. (Order by PB number) . 
Reprints from outside journals and the NBS Journal of Research 
may often be obtained directly from the authors. 
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